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The HIPASS catalogue: Qhi and environmental effects on the H I 
mass function of galaxies 
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ABSTRACT 

We use the catalogue of 4315 extragalactic H I 21-cm emission line detections from HIPASS 
to calculate the most accurate measurement of the Hi mass function (HIMF) of galaxies 
to date. The completeness of the HIPASS sample is well characterised, which enables an 
accurate calculation of space densities. The HIMF is fitted with a Schechter function with 
parameters: a = -1.37 ± 0.03 ± 0.05, log(M*,/M ) = 9.80 ± 0.03 ± 0.03 h^, and 
9* = (6.0 ± 0.8 ± 0.6) x 10~ 3 h% 5 Mpc" 3 (random and systematic uncertainties at 68% CL), 
in good agreement with calculations based on the HIPASS Bright Galaxy Catalogue, which is 
a complete, but smaller, sub-sample of galaxies. The cosmological mass density of H I in the 
local universe is found to be fl m = (3.5 ± 0.4 ± 0.4) x lO - */^ 1 . This large homo geneous 
sample allows us to test whether the shape of the HIMF depends on local galaxy density. We 
find tentative evidence for environmental effects in the sense that the HIMF becomes steeper 
toward higher density regions, ranging from a « —1.2 in the lowest density environments 
to a w —1.5 in the highest density environments probed by this blind H I survey. This effect 
appears stronger when densities are measured on larger scales. 
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1 INTRODUCTION 

The concept of the H I ma ss function (HIMF) of galaxies was in- 
troduced by Brieas 1 1990) as a diagnostic tool for assessing the 
completeness of optical galaxy catalogues against 21-cm Hi sur- 
veys. In that paper the then available luminosity functions and 
an adopted relation between gas richness Mui/L and luminos- 
ity L were used to make the first preliminary HIMF calcula- 
tions. Since then much effort has been put into improving the 
calculations, mostly based on blind 21-cm Hi surveys because 
these - by defini tion - have no bia s against optically faint, gas- 
bearing galaxies jZwaan_eU^ll997|; |Knum^Korte\vea_eL^[l99^ 



Schneider, SuitzalL ^^rjser!bera ~ ^9^nR^senb^rg^^chneMen 
2002; Zwaan et al. 2003, hereafter Z03). The accuracy with which 
the HI MF ca n be evaluated has improved significantly since 
Briggs' ll990l) original calculation, but interestingly, the shape and 
normalisation of the field galaxy HIMF have not changed much. 
This is illustrative of the fact that the galaxy population seen in 
21-cm emission is essentially the same to that seen in the opti- 
cal or IR wavelength, but obviously weighted toward the more gas 
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rich galaxies, and hence to the late morphological types. Two in- 
teresting conclusions can be drawn from the results of the 21-cm 
surveys: 1) There does not seem to be a large population of low 
surface brightness galaxies that is overlooked in optical surveys 
(cf.. lDisnevlll97d) . This conclusion is also reached by deep large 
area optical survey s such as the Sloan Digital Sky Survey (SDSS, 
Blanton et al. 2001 ) and the Two degree Field Galaxy Redshift Sur- 
vey (2dFGRS C ross & Driver 2002), as well as by the Millennium 
Galaxy Survey (MGC, Driver et al. 2005). 2) There appear to be no 
- or at least very few - isolated self-gravitating H I clouds without 
stars, or "dark galaxies". These results are very important for our 
understanding of the galaxy population of the local universe. 



Once an accurate HIMF has been measured, it can be used to 
evaluate the total mass density of neutral hydrogen atoms in the lo- 
cal Universe, Qm(z = 0). At higher redshifts, the highest column 
density QSO absorption line systems are used to measure Qm, and 
these measurements have shown that the comoving neutral gas den- 
sity at z = 2 — 5 is approximately a tenth of a percent of the present 
critical density of the Universe and about h alf of the mass d ensity 
in stars at the pres ent epoch (e.g., Storrie-Lombardi & Wolfe 2000; 
IPeroux et alJl2003T) . For a good understanding of the evolution of 



2 M.A. Zwaan et al. 



the Universe's Hi content, the measurement of a local benchmark 
is obviously very important. 

One aspect of the HIMF, its dependence on environment, has 
not been studied in detail because samples of H I-selected galax- 
ies were not large enough to address this issue. Recently, there 
has been much inter est in the environmental effects on star for- 
mation rate iLewis et al j|2002j : [Baloeh et alJl2004T) . the luminos- 
ity function l Croton et al. 2004), and luminosity, surface bright- 
ness and colour iHo ge et alj 120031) . Since the Hi properties of 
galaxies are intimately tied to the process of star formation, 
gaining knowledge on the environmental effects on Hi proper- 
ties will help understanding the above mentioned results. Two 
previous studies discus sed en vironmental effects on the HIMF: 
lRosenbe~& Schneid er] j2002j) compared the field HIMF w ith that 
in a cluster region and lSpringob. Havnes & Giovanellil 120041) stud- 
ied the HIMF of optically selected galaxies. In this Letter we will 
present for the first time a study of the environmental effects on the 
HIMF for a large sample of H I-selected galaxies. 

The measurement of an accurate HIMF was one of the pri- 
mary science drivers for the H I Parkes All Sky Survey (HIPASS). 
This survey, which blindly searched the whole southern sky for 21- 
cm H I line emission, has now been completed and the catalogue 
of 4315 extragalactic Hi signals has been released to the public 
iMever et all2004tlZwaan et alj2004l) . A HIMF b ased on the 1000 
HIPA SS galaxies with the highest H I peak fluxes l Koribal ski et alJ 
120041) . was published in Z03. This Letter presents the state-of- 
the-art HIMF using the the full HIPASS catalogue. The velocity 
coverage of HIPASS is -1200 to 1 2, 700 km s" 1 and t he typical 
noise level is 13 mjy beam -1 . See lMever et alj 120041) for more 
details on the HIPASS catalogue. Throughout this Letter we use 
Ho = 75kms -1 Mpc -1 to calculate distance dependent quanti- 
ties. 



2 CALCULATING THE HIMF 

In Z03 we derived a bivariate stepwise maximum likelihood 
(2DSWML) technique, which solves for the space density of ob- 
jects as a function of Hi mass (Mm) and velocity width (W) si- 
multaneously. For bins j in log A/hi and bins k in log W the space 
density 6jk can be determined iterat ively by eva luating th e max- 
imum likelihood solution derived by Efstathiou, Ellis, & Peterson 
(1988): 




(1) 



where rijk is the number of galaxies in each bin, Hijk is a function 
that causes the summation to only go over over the galaxies that are 
within the volume accessible to sources in bin jk. The denominator 
Gi is the integral over the bivariate H I mass function defined as 



e, = 



C{S',S' in t)d(Mm,W)dM m dW, (2) 



where C is a function that makes the integral go over the mass and 
velocity range that source i could have and still be detectable at its 
distance Di, and 



S' = S^Mm/M^iWyW), 
SLt = SLM m /M m 



are the effective peak flux and integrated flux of source i if its H I 
mass would change from M m to Mm and its velocity width would 
change from W l to W . 

In Z03 we specifically designed this method for a peak-flux- 
limited sample. However, the equations given above are generally 
applicable for samples constructed using different selection crite- 
ria. The function C, which can be interpreted as a completeness 
function, is defined by the selection criteria of the sample. For ex- 
ample, by setting C = for S' smaller than the peak-flux-limit 
and C = 1 otherwise, the method is again equivalent to the peak- 
flux-limited case in Z03. 

It can be easily seen that the part between the parentheses in 
equationQis essentially the effective volume available to sources 
in bin jk. The important modification that we make here to the 
SWML method, is that instead of evaluating the summation per bin 
jk, we calculate the effective volume for each galaxy individually 
by iteration: 



Y^C(S",S!Ld/@i, 



where S" and S" nt are defined as 

S" = S j {M m /M m ){W j /W l ) 
S'L = S( nt Mm/M m . 



(5) 



(6) 
(7) 



(3) 
(4) 



Stable solutions for V^g are found after ~ 25 iterations. These val- 
ues of VJfj are found for each galaxy individually and are the max- 
imum likelihood equivalents of V ma x for the standard El/V ma x 
method l Schmidt 1968). Summing the values of 1/V^ ff in bins of 
Hi mass gives the original 2DSWML solution for the Hi mass 
function. However, for maximum likelihood estimators of the space 
density of objects, the normalisation is lost and has to be deter- 
mined afterwards. This is done by calculating the mean galaxy den- 
sity n via the minimum variance estimator iDavis & Huchri 19821) 
and setting n equal to the integral over the H I mass function. To 
take into account the effects of sample variance on galaxy den- 
sity, the selection function is weighted by the inverse of the sec- 
ond moment of the two-point correlation function, which is set to 
J 3 = 8000 /if 5 3 Mpc 3 iMever et all2005l) . The selection function 
S is calculated similarly as in Z03. Only the distance range corre- 
sponding to 0.001 < S < 0.1 is used in the calculation of n, to 
avoid adding noise from the farthest and nearest regions of velocity 
space where the number of galaxies is low. After this normalisation 
is applied to , these values represent true volumes. The method 
was tested extensively using synthetic samples, similar to what is 
described in Z03. 



3 RESULTS 

In principle, the 2DSWML technique can be applied to the full 
HICAT sample, including the detections for which the expected 
completeness is very low. However, including the detections for 
which C « 1 implies large corrections to the expected space 
densities (see Eq|5J. To avoid this, we choose to cut the sample at 
C > 0.5. This reduces the total sample by only 7% to 4010. For 
the completeness function C we adopt the parametric approxima- 
tion based on an error function calculated in lZwaan et alj 120041) . 
and to account for the varying noise level across the sky, we adopt 
the completeness correction suggested in that paper. We note that 
the latter correctio n has an insignificant effect on the HIMF. In 
IZwaan et all <2004 we found that HICAT is complete at the 95% 



HIPASS H I mass function 3 




7 8 9 10 11 

og M H! [M s ] 



Figure 1. Top: H I mass function from HICAT. The solid line is a Schechter 
fit to the points, the best-fit parameters are shown in the lower left corner; 
Bottom: Distribution of H I masses used for the HIMF calculation. 

(99%) level at a peak flux of 68 (84) mJy and at an integrated flux 
of 7.4 (9.4) Jykms" 1 

FigureQshows the result of the analysis. The top panel shows 
the HIMF, the bottom panel shows a histogram of the number of 
galaxies in each H I mass bin. The errors on the HIMF points are 
from Poisson statistics. A Schechter function defined by 

e(« H ,M«», = «-(^)"«p(-^),(^), 

is fitted to the points and shown by a solid line. The parameters of 
this fit are: a = -1.37±0.03, log^^/M©) = 9.80±0.03 h^, 
and 6* = (6.0 ± 0.8) x 10~ 3 h% Mpc~ 3 . The uncertainties on 
the Schechter parameters are determined by jackknife resampling 
jLuptonll993l) . where the sky is divided into 24 equal area regions. 
This jackknife resampling procedure includes the uncertainties due 
to large-scale structure and possible errors due to the varying noise 
level across the sky. The dimensionless covariance between a and 
log A/hi is r(a, log Mhi) = —0.53, indicating the strong correla- 
tion between these two parameters, as is normally seen in galaxy 
luminosity and mass function fits. We note that log Mj^ and 8* are 
also highly covariant: r (log A/jjr, 6*) = —0.83, which illustrates 
that quoting values of log A/j^i or 9* individually is not very mean- 
ingful: the combination of these two parameters defines Qui- For 
completeness, we note that r(a, 8*) = 0.32. 

For comparison, the HIMF derived from the BGC (Z03) is 
plotted as a dashed curve. The new HICAT HIMF is slightly steeper 
than the Z03 HIMF, for which we found a = -1.30 ± 0.08, 
but within the errors. The new measurement is also in reasonable 
agreement with previous HIMF measurement s from blind 21 -cm 
surveys, which showed values of a = —1.2 (Zwaanetal. 1997), 
a = -1.52 jKilborrl200fl. a = -1.51 iHenning et all2000h . and 
q = —1.53 (Rosenbere & Schneider 2002). We refer the reader 
to Z03 for a complete discussion of previous HIMF measurements 
from blind 21 -cm surveys. 

In Figure [2] we show the redshift distribution of the HICAT 
sample. A good method to test the validity of the method used to 
calculate the HIMF is to compare this distribution with the pre- 
dicted distribution from the derived selection function S. This pre- 
dicted distribution is shown by the dashed line, which appears to be 
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Figure 2. Top: H I masses of HICAT detections as a function of their red- 
shift. Points above the solid line are integrated flux-limited at S; nt = 
9.4 Jy km s — 1 <Zwaan et all2004l) . Only points that are used in the HIMF 
calculation (C > 0.5) are shown. Bottom: Redshift distribution of HICAT 
galaxies. The histogram shows the measured distribution and the dashed 
line is predicted distribution based on the selection function calculated from 
the maximum likelihood method. 

an excellent approximation to the real HICAT selection. Similarly 
to what we found in Z03, the number counts predicted from the 
selection function (0.174 deg -2 ) is only 7% lower than the mea- 
sured number counts (0.186 deg -2 ). This difference is due to the 
two over-densities seen at ~ 1500 and ~ 2500 kms -1 . 

3.1 Biases in the HIMF determination 

In Z03 we discussed in detail the various biases in the determi- 
nation of the HIMF based on the Bright Galaxy Catalogue. It 
was found that the selection bias and H I self-absorption have the 
strongest effects on the HIMF and Q,m- For brevity, we choose 
not to repeat here all the tests of biases in the determination of the 
HIMF, but only briefly discuss the ones that may be particularly 
relevant to HICAT and refer the reader to Z03 for a detailed de- 
scription of biases. 

Selection bias — The selection bias is the effect that noise 
on detection spectra preferentially causes peak fluxes to be over- 
estimated, which may result in an increased number of low-flux 
detections in the sample. Similarly to Z03, we simulate this effect 
and find that due to the selection bias the global space density of 
objects is probably overestimated by ~ 10%, and the slope of the 
HIMF becomes steeper by ~ —0.05. Together, this causes an over- 
estimation of Q,m of ~15%. 

H I self absorption — The effect of H I self-absorption is very 
difficult to ass ess from the HIPASS data. We conservatively adopt 
the result from lZwaan et alJ 1 19971) that H I self-absorption causes 
a small decline in the measured space density and the knee of the 
HIMF, resulting in underestimation of the measured flm of ^15%. 

Deviations from Hu bble flow — 

iMasters. Havnes. & Giovanellil 120041) recently argued that 
calculating distances to galaxies by assuming a pure Hubble flow 
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instead taking into account the local velocity field, can lead to 
errors in the slope of the HIMF. We test this by adopting the 
Moul d*etal] l2000l) distance model and recalculate the HIMF. 
The low-mass slope of this HIMF is a — —1.39 ± 0.03 not 
significantly different from our HIMF based on pure Hubbl e flow 
distances in the Local Group reference frame. [Masters et all 1 2004) 
concluded that the Z03 HIMF slope was probably underestimated 
and they suggest a slope of a = —1.4 ± 0.1, which is in excellent 
agreement with our new HICAT slope. The median distance of the 
HICAT sample is larger than that of the Z03 sample, which implies 
that the fractional distance errors due to peculiar motions are 
smaller. This explains why taking into account the local peculiar 
velocity field does not have a significant effect on the HICAT 
HIMF slope. 

Lower distance limit — To avoid confusion with Galactic 
high-velocity clouds (HVCs), we applied a lower limit of v gsI < 
300 km s -1 to HICAT. Obviously, this constraint also excludes 
very nearby, low mass galaxies from the present analysis. Our pre- 
vious HIMF determination based on the Bright Galaxy Catalogue 
(Z03) had a measurement at M H i = 1O 6 ' 8 M , but the lowest Hi 
mass point for HICAT is at Mm = 1O 7 ' 2 M . In Z03 we argued 
that for peak-flux-limited samples, imposing a lower distance limit 
can cause a slight drop of the HIMF at very low Hi masses. For 
HICAT we find that this effect is not measurable because the sam- 
ple is not strictly peak-flux-limited. 



4 ENVIRONMENTAL EFFECTS ON THE HIMF 

lRosenbe~ & Schneider (2002) speculated on the effects of envi- 
ronment on the slope of the HIMF. They concluded that the HIMF 
in the immediate region around the Virgo cluster was a f» —1.2, 
flatter than in the field region where they found a ~ —1.5. Our 
sample of HICAT galaxies is sufficiently large to test the effect of 
local density on the shape of the HIMF, by subdividing the cata- 
logue into samples drawn from different local densities. 

We choose to parameterise local density by calculating the 
distance r n to the n-lh nearest neighbour for all HICAT galaxies. 
Densities are calculated as p n = n/[(4/3)7rr^], and normalised to 
the median density to obtain p„/ < p„ >. We divide the HICAT 
sample into five logarithmically spaced subsamples with different 
mean densities. Since HICAT is increasingly sparsely sampled at 
larger redshift, nearest neighbour measurements increase rapidly as 
a function of redshift. Therefore, different subsamples could effec- 
tively be drawn from different redshift ranges. This, in turn, implies 
preferential selection of certain H I mass ranges for each subsam- 
ple (see Fig.|2j, which could bias the HIMF calculation. We cor- 
rect for this effect by first fitting a power-law to the conditional 
probability distribution of nearest neighbour distance as a function 
of redshift. Next we subtract this power-law fit from the nearest 
neighbour distances to obtain an unbiased nearest neighbour mea- 
surement of which, however, the absolute calibration is now lost. 
We test that our corrected nearest neighbour measurements sample 
all distances equally by calculating for each subsample the mean 
and the dispersion of the redshifts. We confirm that there are no 
significant differences in the redshift distributions of the subsam- 
ples. Finally, to prevent including galaxies from the most sparsely 
sampled redshift range, we limit our calculations to galaxies with 
v gBI < 5000kms~\ 

The HIMF as a function of local density ps/ < ps > is shown 
in Fig. [3] We find evidence for a weak dependence of the HIMF 
slope on density environment, in the sense that the low-mass end of 
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Figure 3. The H I mass function as a function of local galaxy density. The 
HICAT sample has been divided into five logarithmically spaced subsam- 
ples according to the density out to the fifth nearest neighbour. The HIMF 
for each subsample is shown with the corresponding best bit Schechter 
function. The HIMFs are offset from each other vertically by 0.5 dex. 
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Figure 4. The low-mass slope a of the Schechter function fits to the H I 
mass function as a function of local galaxy density. The different panels 
show the results based on the first, third, fifth, and tenth nearest neighbour, 
from left to right. 

the HIMF becomes steeper in higher density regions. There is no 
corresponding correlation between the density and the characteris- 
tic mass Afj^j. In addition to this analysis based on the fifth nearest 
neighbour, we also calculated densities on smaller and larger scales, 
out to the first, third, and tenth nearest neighbour, respectively. The 
results based on these measurements are shown in Figure|4| where 
the Schechter low-mass slope a is plotted as a function of density. 
Interestingly, the effect of steeper HIMFs toward higher density re- 
gions appears to become stronger toward larger scales, indicating 
that the environmental effects on the H I properties of galaxies are 
not restricted to short distance effects such as galaxy interactions, 
but also operate on larger scales. 

Our observed trend between the slope of HIMF l ow-mass en d 
with density is in qualitative agreement with lCroton et alJ 120041) . 
who use the 2dF galaxy redshift survey to show that the slope of 
the luminosity function is weakly dependent on envi ronment. Our 
results cannot be quantitatively compared to those of ICroton et alJ 
120041) because we use a different method to parametrise local 
density. Furthermore, our de nsity measurements are based on H I- 
selected galaxies only, while Croton et al. 1 2004) use optically se- 
lected galaxies. It should also be kept in mind that the dynamic 
range in densities probed by HIPASS is sm aller th an that probed by 
optical surveys. For example, Waueh et alJl2002h showed that clus- 
ter centres that are readily apparent in the counts of optical galax- 
ies are only detected as over-densities of a few times the density of 
H I-detected galaxies in the field. At first sight, our findings seem to 
be in conflict with the recent results of Springob et al. (2004), who, 
based on optically selected galaxies, claim a flattening of the HIMF 
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toward higher density regions. However, the statistical significance 
of their result is low and the HIMFs from the highest and lowest 
density subsamples are consistent at the la level. 

Interestingly, recent 2dF and SDSS results show an anti- 
correlation between star formation rates based on the equivalent 
width of the Ha emission line and galaxy density llxwis et alj 
l2002tlG6mez et alj2003HBalogh et alj2()04l). o ut to l arge distances 
from cluster cores. [Kauffmann et alj 120041) and iTanaka et alJ 
120041) also find that the specific star formation rate SFR/M,, is 
strongly dependent on environment, but furthermore conclude that 
this relation is most pronounced for the lower mass galaxies. This 
implies that the comparative rate at which lower mass galaxies 
consume their innate Hi gas with respect to higher mass galax- 
ies is greater in low density environments than high density envi- 
ronments. Hence, this effect leads to a flattening of the HIMF in 
low density regions. Of course, this is a simplified view of galaxy 
evolution because it is well established that gas accretion plays 
an important role in the tot al neutral gas budget of galaxies (e.g., 
Matte ucci & FrancoiJl989() . but it can go a long way in explaining 
our observed correlation between HIMF slope and density. 

An alternative explanation can be found in the halo- 
occupation model in which galaxy properties depend solely on 
dark matter halo mass, which is predicted by the cold dark mat- 
ter (CDM) theory. iMo et alJ <2004 predict that in this model the 
shape of the luminosity function of late type galaxies is indepen- 
dent of local density because these galaxies live in low-mass halos. 
Early type galaxies mainly exist in high-mass halos, which are as- 
sociated with higher density regions. It is known that the H I masses 
of early type galaxies are typi cally low, and thei r HIMF is steeper 
than that of late type galaxies IZwaan et alj2 003). Therefore, a rel- 
atively higher contribution of lower H I masses (and hence a steeper 
HIMF) in higher density regions is a direct prediction of the CDM 
model. 



5 CONCLUSIONS 

In Z03 we presented the Hi mass function of galaxies based on 
the 1000 HIPASS galaxies with the highest peak fluxes. That pa- 
per gives a detailed description of the HIMF calculation, the pos- 
sible effects that may bias the calculation, and discussed the re- 
sults. This Letter presents an updated HIMF measurement from 
the full HIPASS sample. We find a = -1.37 ± 0.03 ± 0.05, 
log(iW"Hi/M©) = 9.80 ± 0.03 ± 0.03 h^, and 6* = (6.0 ± 
0.8±0.6) x 10 -3 /175 Mpc -3 (random and systematic uncertainties 
at 68% confidence), in good agreement with the Z03 results. The 
slope of the HIMF is comparable to that of the optical luminosity 
function from the SDSS measured down to low absolute magni- 
tudes of M r ~ -13 jBlanton et all2004 . 

We find tentative evidence for a steepening of the low-mass 
end of the HIMF as a function of local galaxy density. This effect 
becomes stronger if local densities are measured over larger scales, 
which indicates that environmental effects on the H I properties of 
galaxies are not restricted to short distance effects such as galaxy 
interactions. 

HIPASS is a relatively shallow survey compared to some of 
the small area deep surveys. However, because of the large sky cov- 
erage the number of detections is very large and the HIMF can be 
measured with high accuracy. HIPASS has particularly good statis- 
tics at the high mass end of the HIMF, which dominates the mass 
density of neutral hydrogen. Therefore, fim is measured with high 
accuracy: Q H l = (3.5 ± 0.4 ± 0.4) x lO -4 ^ 1 . 



The faint-end tail of the HIMF is measured well down to H I 
masses of 1O 7 ' 2 M0. No other survey has been able to constrain 
the HIMF better at this mass limit. The space density of objects 
much below A/hi = 10 7 Mq is still not accurately known. A num- 
ber of surveys have been able to put upper limits on the space 
densit y of these low Hi mass sy stems in higher density regions 
(e.g Jde Blok et all2002tlZwaanl 2001), but large scale 'blind' sur- 
veys with existing technology are probably not able to make sig- 
nificantly better measurements in this region of parameter space. 
Future large radio telescopes are required to measure reliably the 
space density of objects much below A/hi < 10 7 M@. 
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